Intrinsic Domain Wall Resistance in Ferromagnetic Semiconductors 
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Transport through zincblende magnetic semiconductors with magnetic domain walls is studied 
theoretically. We show that these magnetic domain walls have an intrinsic resistance due to the 
effective hole spin-orbit interaction. The intrinsic resistance is independent of the domain wall 
shape and width, and survives the adiabatic limit. For typical parameters, the intrinsic domain wall 
resistance is comparable to the Sharvin resistance and should be experimentally measurable. 



Understanding magnetic topological defects is crucial 
in developing devices that utilize the electron spin. Do- 
main walls are topological defects between homogeneous 
magnetic domains. The domain wall dynamics has tradi- 
tionally been induced by external magnetic fields. There 
has recently been a large interest in the scientific commu- 
nity on current induced magnetization dynamics, where 
domain walls and domains change in response to applied 
electric currents 0. Domain walls can also be electri- 
cally detected, by their electrical resistance. Knowledge 
of the domain wall's effect on the electrical resistance is 
important for the understanding of spin transport in con- 
densed matter and for the electrical detection of magnetic 
topological defects. 

Transport through domain walls have been extensively 
studied in metallic systems, theoretically and exper- 
imentally The domain wall resistance is defined as 
R w = R — Ro, where R and Rq are resistances with a 
domain wall and with homogeneous magnetization, re- 
spectively. When the domain wall is thinner than the 
mean free path, in the ballistic regime, R w is positive. 
In diffusive systems, when the domain wall is wider than 
the mean free path, the sign of the domain wall resistance 
is still under debate, i.e. the domain wall can increase 
or decrease the resistance of the ferromagnet. In ballistic 
and diffusive metal systems, R w approaches zero with in- 
creasing domain wall width and vanishes in the adiabatic 
limit when the domain wall is much wider than the Fermi 
wavelength. 

Ferromagnetic semiconductors integrate magnetiza- 
tion controlled spin transport with gate controlled carrier 
densities in semiconductors. Domain walls in these sys- 
tems have been recently studied experimentally 0, 0, El 
7] . The strong interaction between the spin of the effec- 
tive holes and their orbits in dilute magnetic semiconduc- 
tors changes the transport properties of magnetic domain 
walls qualitatively. In this Letter, we show that domain 
walls in zincblende magnetic semiconductors have an in- 
trinsic resistance R^ which is the part of R w that sur- 
vives the adiabatic limit. R^ is independent of the width 
and detailed shape of the domain walls and is due to the 
effective spin-orbit interaction. 

Related manifestations of the coupling between the 
spin-orbit interactions and the magnetizations are the 
anisotropic magneto resistance (AMR) H,EJE1EI and 



the tunneling anisotropic magneto resistance (TAMR) 
[jllTi! fl3|. In domain walls, some carriers are prevented 
by the spin-orbit interaction to adiabatically adapt to the 
change in the direction of the local magnetization produc- 
ing a finite R w even in the adiabatic limit. R^ depends 
mainly on three material parameters: the effective spin- 
orbit coupling, the exchange field and the Fermi energy. 
Hole transport in dilute magnetic semiconductors with 
strong spin-orbit interaction is usually difficult to treat 
analytically. Nevertheless, we find an em analytical so- 
lution for the adiabatic domain wall conductance. 

We model hole transport in zincblende magnetic semi- 
conductors by the following dimcnsionless Hamiltonian 



H = aiPiPi - a 2 (JiJjP l Pj) + Jihi 



(1) 



where the subscripts (i,j) = x,y,z and the Einstein sum 
convention is assumed. Furthermore, h(r) is the di- 
mcnsionless exchange field, describing the interaction be- 
tween holes and localized magnetic moments, |h(r)| = 
ho. In Eq. ^ p is the dimensionless momentum op- 
erator and J denotes the angular momentum operator 
for J = 3/2 spins. The parameters cti and a 2 are 
controlled by the hole effective mass and the strength 
of the spin-orbit interaction, respectively. Their rela- 
tion to the Luttinger parameters |l5| 71, 72 and 73 are 
ai = (71 + 572/2)7(71-272) and a 2 = 272/(71-272). 
In our notation, energies, momentums and lengths are 
measured in units of the Fermi energy, Fermi momen- 
tum and Fermi wavelength of heavy holes for a given 
doping in the original non-magnetic host system, /io = 0. 
A six band version of Eq. ^ explains man y fe atures of 
zincblende magnetic semiconducto rs I lfil ll7L . We em- 
ploy the spherical approximation yjj and disregard two 
split-off bands. Thus, 73 =72 and the Fermi energy Ep 
is assumed to be smaller than the split-off energy. We 
compute the conductance G = l/R from the Landauer- 
Biittiker formula G = (e 2 /h)J2nm \t n m\ 2 , where t nm is 
the unit flux normalized transmission between the n and 
m transverse waveguide modes. 

Let us first prove that holes governed by the Hamil- 
tonian Q exhibit an intrinsic adiabatic domain wall re- 
sistance. We separate the Hamiltonian into an intrinsic 
contribution Hi and a collision contribution He by trans- 
forming H into a local basis where the quantization axis 
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for J is parallel to h(r), H ^UHU' 1 =H I + H C - 

Hi = uiPiPi - ct 2 (J l J j p i p ) + Jzh , (2) 
H c = (otiSij - a 2 JiJj)[K i K :j + Ay + 2mpj] , 

where U(r) is a 4 x 4 unitary operator defined such that 
U[h(r) ■ jjc/- 1 = J z h Q . Furthermore, J* = UJiU' 1 , 
Ki = UpiU and Ay = pJJpjU^ 1 are 4x4 matrices 
that do not operate on the spatial coordinate r. 
Consider the adiabatic limit when the width of the do- 
main wall X w is much larger than the Fermi wavelength 
Xp. Hence, Ki ~ 1/X W — ► 0, Ay ~ V^ui ~^ and thus 
He — * 0. Without the effective spin-orbit coupling, 
o;2 = 0, flj is independent of directional variations in 
the exchange field and the domain wall resistance van- 
ishes, R w = 0. In contrast, with finite effective spin-orbit 
coupling, a 2 ^ 0, Hi varies for systems with or with- 
out variations in the exchange field, since Ji differs from 
Jj. In other words, the spin polarization of the transport 
channels is anisotropic, giving rise to a finite intrinsic 
domain wall resistance. 

In order to quantify our findings, we consider the lin- 
ear response of a rectangular conductor with a cross 
section A — L X L Z and periodic boundary condition 
(PBC) along the x and z directions. Transport is along 
the y-axis. We consider three types of domain walls: 
Bloch ZX-wall, Neel ZY wall and Neel YZ-wall de- 
scribed by Hy) = [fc(y),0,f 1 (y)], h(y) = [0,f 2 (y),hm 
and h(y) = [0, f 2 {y)], respectively. Here, fi(y) = 

h tanh(y/X w ) and f 2 (y) = h / cosh(y/X w ). 

Translation invariance conserves the transverse mo- 
menta k x and k z that label the transport channels. Ad- 
ditionally, each (k x ,k z ) channel contains four internal 
spin channels originating from the four spin-orbit cou- 
pled bands. The distribution of the transport channels 
at position y, T y {k x ,k z ), may be found by solving the 
eigenvalue problem H(k x ,y,k z )ip y = Eptp y . The Fermi 
energy, Ep, is position independent. Due to the inter- 
play of the spin-orbit interaction and the magnetization, 
the number of open channels is anisotropic in the mo- 
mentum k x -k z space, unlike systems where the effective 
spin-orbit interaction vanishes. In the adiabatic limit, a 
conducting channel must exists throughout the system in 
order to contribute to the conductance. Mathematically 
this can be expressed as G — {e 2 /h) J2k x T(k x , k y ), 
where T{k x ,k y ) — C\ y < L_ OQ T y (k x ,k z ) is the intersection 
of all distributions of transport channels as one traverses 
the domain wall, see Fig.l. In the case of a homoge- 
neous ferromagnet, the conductance is Go = \/Rq = 
(e 2 /^) Efc^fcz Ty=-oa(k Xl k z ). If and only if T y= -oo is 
identical to or is a subset of T y for all y then = 0. 
This occurs only when the effective spin-orbit interaction 
or ho vanishes. 

For Bloch walls, T y (k x ,k z ) is conformal under trans- 
lation and simply rotates in the same direction as the 
exchange field, see Fig.l. Traversing from y — -co to 



oo corresponds to a 7r-rotation of Ty = -oo(k x , k z ) forcing 
T(k x , k y ) to be circular symmetric, see Fig.l. An analyt- 
ical expression for the ballistic conductance through an 
adiabatic Bloch wall can be written as 

° B = i Ke ( fcl ) 2 + Ke ( fc 2) 2 + ^ e ( fc 3) 2 + 5Re(fc 4 ) 2 }(3) 

where 5fte denotes the real part and 

fci = ^J{2 ll Ep + D)B , fc 4 = V(27i^f - D)B 

k 2 = H(Xi)V^f-3/i /2 

+ H(-Xi) V(7i ~ 2 72 )(2£ F + h )B 

k 3 = E{X 1 )^/( 11 -2 l2 )(2E F + h )B 

+ S(-Xi) Z(X 2 )^Ep~ 3h /2 

+ H(-X 2 )v/(7i " 2j 2 )(2E F - h )B 

where is the Heaviside step function with 2(0) = 

0.5, X l = (2 72 ^ i r)/(7 1 +7 2 ) ~ and X 2 = (4 72 £ J r)/( 7l + 
472) — h . Furthermore, B= 1/(27! + 4 72 ) and 

D = /i (7i ~ 272) + 

2^0(7? ~ 7i72 ~ 2 7 |) - 2h 0l2 E F ( 11 - 72 ) + 4r,*E* . 

For Neel walls, T y (k x ,k z ) is not conformal under trans- 
lation and Eq. |31 does not apply. 

Beyond the adiabatic app roximation, we use a stable 
transfer matrix method |20| to find the conductance nu- 
merically. The system is discretized into N x x N y x N z 
lattice points with lattice constants a x , a y and a z . We use 
L Z = L X = N x a x = 5, a x =a z = 0.05, L y = 150, a y = 0.08 
and 7! = 6.8. Varying 72 or ho changes the reservoir 
properties. We have chosen to vary the Fermi energy, 
Ep{^ 2 ,ho), keeping the hole carrier density constant, to 
mimic experimental conditions. Ep("f 2 , ho) is calculated 
numerically. We fixed the volume density of holes such 
that E F (j 2 = 2.1, ho = 0) = 1. Experimentally [U Ej, 
Ep(h )/h ~ 1 which corresponds to ho ~ 1.5. The con- 
vergence of the numerical scheme is verified by check- 
ing the unitarity of the scattering matrix. We present 
our results in terms of relative domain wall resistance, 
$tw = (R — R )/R Q . 

For increasing A^, the direct collision between holes 
and the domain wall gradient decreases and %t w (X w ) 
drops rapidly to its intrinsic value, Fig.l. Note that 
^w(X w ) is very close to its intrinsic value already when 
X w ~ 1 for ZX and ZY walls, and X w ^0.5 for a YZ wall 
e.g. when the domain wall width is equal or smaller than 
the Fermi wavelength of heavy holes. 

The relative intrinsic domain wall resistance as a func- 
tion of the exchange field is shown in Fig. 2. Generally, 
Sffc^/jo) increases for increasing ho, due to an increase in 
the anisotropy in T y (k x ,k z ). For the YZ wall, however, 
the Sharvin resistance i?o increases even faster than R 
leading to a weak reduction in at large ho- We also 
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FIG. 1: Relative domain wall resistance versus wall width for 
ZX wall (triangle), ZY wall (circle), and YZ-wall (square). 
The solid line displays the intrinsic resistance derived from 
Eq. |3 Bottom right inset:, $l w for a film shaped system 
with Dirichlet boundary conditions. Top left to right insets: 
Distribution of conducting channels T y (k x ,k z ) for y = — oo, 
y = —0.88\ w , y~0, and the intersection T(k x ,k z ), for a ZX 
wall. 




FIG. 2: Intrinsic relative domain wall resistance versus ex- 
change field for ZX-wall (triangle), ZY-wall (circle), and YZ- 
wall (square). The solid line displays the analytical result for 
a ZX-wall, Eq. (3). Inset: 5R^,(/io) f° r a film shaped system 
with Dirichlet boundary conditions. 



see in Fig. 2 that the numerical result for the Bloch wall is, 
within 2%, identical to the analytical result derived from 
Eq. |3 for ho > 0.5. There are small deviations at small 
ho < 0.5. Here, the adiabatic limit requires very large 
A-U, and system sizes larger than one we used, L y = 150, 
for convergence. 

The relative intrinsic domain wall resistance as a func- 




FIG. 3: Intrinsic relative domain- wall resistance as functions 
of hole spin-orbit coupling parameter, ^(72): ZX-wall (tri- 
angle), ZY-wall (circle) and YZ-wall (square). Solid line 
shows the analytical result for a ZX-wall. Inset: 3^(72) for a 
film shaped system with Dirichlet boundary conditions. 



tion of the effective spin-orbit coupling is shown in Fig. 
3. 5i4(72) increases monotonically with increasing 72 for 
all walls. This is due to an increase in the anisotropy 
in T y (k x ,k z ). We also see that the numerical result for 
the Bloch wall is identical to the analytical result de- 
rived from Eq. [3] for 72 < 2.7. The small deviation at 
large 72 > 2.7 is caused by the same effect as previously 
discussed for Jf^C^O < 0.5). Here, increasing 72 reduces 
the effective amplitude of the exchange field along certain 
directions. 

For reasonable parameters 72 ~ 2.1 and ho ~ 1.5, 
5R^~0.9 and 0.7 for the ZX and ZY walls, respectively. 
In other words, a domain wall removes nearly half of the 
Sharvin conducting channels in ballistic adiabatic trans- 
port, an effect that should be clearly measurable. An 
interesting question is how much of the intrinsic resis- 
tance still remains in the diffusive transport regime. We 
know that the anisotropy in the distribution of conduct- 
ing channels with respect to the direction of the exchange 
field do survive in the diffusive limit leading to the AMR 
effect 0,111, El EI- We therefore expect that at least part 
of the intrinsic resistance will also survive the diffusive 
limit. Furthermore, the intrinsic domain wall resistance 
is reached already for domain wall widths comparable to 
the Fermi wavelength. Consequently, we expect that the 
ballistic intrinsic domain wall resistance will be impor- 
tant even in rather dirty state-of-the art dilute ferromag- 
netic semiconductors where the mean free paths are not 
very much smaller than the Fermi wavelengths. 

Let us explain the rigidity of the intrinsic domain wall 
resistance against variations in the domain wall's width 
and shape using topological arguments. Consider a Bloch 
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ZX wall. First, define the order parameter space as a 
two dimensional space in which two dimensional vectors 
map to points 22], e.g. the real space vector (h x x+h z z) 
maps into the point (M X ,M Z ) — (h x ,h z ) in the order 
parameter space. Second, note that the distribution of 
conducting channels, T y (k x ,k z ), depends on the average 
exchange field within a wave packet h(y) — + h(y). Thus, 
the conductance G ~ C\ y < L_ 00 T y (k x , k z ) depends on h(y) 
or, more precisely, its mapping in the order parameter 
space. For X w <C Xp, h(y) maps, approximately, into a 
half ellipse, M 2 X 2 F /9X 2 W +M 2 = h 2 where M x > 0. The 
mapped curve becomes more circular for increasing A^,. 
For A™ > Xp, h(y) maps to a half circle, M 2j rM 2 — h\ 
where M x > 0. In this smooth (adiabatic) limit, ZX 
domain walls of all widths and shapes map into the same 
half circle, and thus have the same conductance. 

So far, we have assumed periodic boundary condition 
in the transverse directions. Spin and orbital motion 
are coupled at the boundaries due to the spin-orbit in- 
teraction, which could change the results. In order to 
address this question, we have also numerically com- 
puted the conductance using Eq.Q]with Dirichlet bound- 
ary condition, ip(0,y, z) = tp(L x + a x ,y, z) = ip(x,y,0) = 
ip(x, y, L z + a z ) = 0. To mimic experimental available 3D 
film shaped conductors, we use L x = 6, a x = 1/3, L z = 3, 
a z = 3/10 and L y = 100, a y — 1/Att. It turns out that 
the Dirichlet boundary condition together with small L z 
prevent the development of extreme anisotropics in the 
distribution of conducting channels, expected at large ho 
and 72. This leads to a reduction of Si^ by a factor of 
2 for the ZX wall. 3^ for the ZY and YZ walls are less 
affected since the distribution of conducting channels for 
h||y is isotropic. The relative domain wall resistance as 
a function of X w is shown in the inset of Fig. ^ Simi- 
lar to the PBC case, 3?™ (A) drops to its intrinsic value 
at A«; ~ 1 for ZX and ZY walls, and Au, ~ 0.5 for YZ 
wall. 3?^(/io) is shown in the inset of Fig. 2. Here, 3?^, 
for the ZX-wall forms a plateau for ho > 0.5 above which 
the anisotropy in the conducting channels is prevented to 
develop further. We see in the inset of Fig.[3]that 3?^ (72) 
for the ZY and ZX walls develop a peak/plateau around 
72 ~ 2 where the anisotropy in the conducting channels 
is prevented to develop further. In any case, the intrinsic 
resistance remains well defined and finite in a 3D film 
shaped conductor with Dirichlet boundary condition. 

In conclusion, we have shown that dilute magnetic 
semiconductors exhibits an intrinsic domain wall resis- 
tance. The intrinsic domain wall resistance depends only 
on the map of the domain wall in the magnetic order 
parameter space and not on their real space width and 
shape. An analytical expression for the adiabatic con- 
ductance for Bloch walls is given. For typical parame- 
ters, the intrinsic domain wall resistance is comparable 
to the Sharvin resistance, and should therefore be exper- 
imentally measurable. The domain wall resistance drops 



to its intrinsic value when X w approaches half the Fermi 
wavelength for heavy holes or a quarter of Xf for light 
holes. These values are not too far from mean free path 
accessible presently accessible, e.g. in Ref. 4]. From 
general topological arguments, we show that topologi- 
cal magnetic defects have an intrinsic resistance against 
transport of carriers with strong spin-orbit coupling that 
survive the adiabatic limit. 
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